Stretching of polymers around the Kolmogorov scale in a turbulent shear flow 
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We present numerical studies of stretching of Hookean dumbbells in a turbulent Navier-Stokes 
flow with a linear mean profile, (u x ) = Sy. In addition to the turbulence features beyond the viscous 
Kolmogorov scale n, the dynamics at the equilibrium extension of the dumbbells significantly below 
rj is well resolved. The variation of the constant shear rate S causes a change of the turbulent velocity 
fluctuations on all scales and thus of the intensity of local stretching rate of the advecting flow. The 
latter is measured by the maximum Lyapunov exponent Ai which is found to increase as Ai ~ S 3 ^ 2 , 
in agreement with a dimensional argument. The ensemble of up to 2 x fO 6 passively advected 
' dumbbells is advanced by Brownian dynamics simulations in combination with a pseudospectral 

integration for the turbulent shear flow. Anisotropy of stretching is quantified by the statistics 
of the azimuthal angle <f) which measures the alignment with the mean flow axis in the x-y shear 
(N ■ plane, and the polar angle 8 which determines the orientation with respect to the shear plane. 

The asymmetry of the probability density function (PDF) of 4> increases with growing shear rate 
5*. Furthermore, the PDF becomes increasingly peaked around mean flow direction {<j> = 0). In 
contrast, the PDF of the polar angle 8 is symmetric and less sensitive to changes of S. 
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I. INTRODUCTION 



When a few parts per million in weight of long-chained polymers are added to a turbulent fluid its properties change 
drastically and a significant reduction of turbulent drag is observed (see e.g. Refs. 1 and 2 for reviews). Although the 
phenomenon is known from pipe flow experiments for almost 60 years Q a complete understanding is still lacking. It 
is observed in wall-bounded flows 0. 0.[f| as well as in bulk turbulence. 0, The macroscopic drag reduction seems 
to go in line with a microscopic transition of the macromolecules from a preferentially coiled near-the-equilibrium to 
a stretched non-equilibrium state. The relation of the polymer dynamics at small scales to the flow statistics at larger 
f~>) \ scales is therefore essential for solving the drag reduction problem. [alSl 

The dynamics of single long polymer chains in simple and partly steady flows is already quite complicated and rich 
t-H | with respect to conformational variety. This was demonstrated in experiments and Brownian dynamics simulations. 

|ToL ITU Il2l ITli. ITU ] Only recently, a more complex situation was investigated - the stretching of polymer chains 
in random isotropic and random linear shear flows - by means of experiments in the so-called elastic turbulence 
limit. pL5L Il6| A naly tic solutions are possible for the kinematic stretching of polymers in isotropic and white-in-time 
Kraichnan flow. 11 71 H E HJ The polymer chain is there described as two beads that are connected by a spring. 
i^H ■ The model is known as the mesoscopic dumbbell model and the entropic spring force follows either a linear Hookean 
Ch ■ or a nonlinear law. The latter captures the finite extensibility of the polymers. |2l| Qualitatively new effects can be 
observed in a shear flow such as the tumbling, a flipping of the chain in the plane that is spanned by the shear 
flow. 0,|23 Analytical predictions for the tumbling statistics of dumbbells are possible when the isotropic Kraichnan 
flow is superposed with a linear shear flow at a very large shear rate. [23 The dynamics of the angular degrees of 
freedom can then be simplified significantly because the large shear aligns the dumbbells preferentially with the mean 
flow direction. 

Although these studies give us a lot of helpful insights about the statistics of the extension and orientation of the 
polymer chains, the situation in a turbulent Navier-Stokes fluid is more complex. The flow structures are correlated 
in space and time and it can be expected that they cause significant differences in the stretching history of a chain 
while moving through the flow. Furthermore, turbulent fluctuations and the large scale shear are intimately coupled 
to each other. In other words, one cannot consider them as independently adjustable parameters. An increase of 
the shear rate S will change the Reynolds number of the flow and more importantly the ratio and magnitude of the 
turbulent velocity fluctuations (uf) with i = x,y,z (see Refs. 24 and 25). For example, the presence of shear causes 
the generation of streamwise streaks, that enlarge the fluctuations in the streamwise direction. [2(1 As a consequence, 
we find that the maximum stretching rate grows with respect to S faster as recently predicted. 23] 



* Present address: International Centre for Theoretical Physics, 34014 Trieste, Italy 

t Present address: Department of Mechanical Engineering, Ilmenau University of Technology, D-98684 Ilmenau, Germany 



2 



The present work consists of two parts. The first part will focus on the statistics of the flow itself, at scales 
below and above the viscous Kolmogorov scale, r\ = aai/4 > with v being the kinematic viscosity and (e) the mean 
energy dissipation rate of the flow. We consider a simple turbulent shear flow with a linear mean profile which 
allows for studying the effects of shear on Lagrangian stretching rates. This in turn will reveal a different regime of 
polymer stretching in comparison with an analytic model of Chertkov et al. |23j] The second part presents studies of 
the kinematics of polymer stretching in such flows. We will demonstrate that the interplay between the shear rate S 
and the Weissenberg number Wi is important for the stretching statistics. The dimensionless parameter Wi relates 
the local stretching rate of the flow to the relaxation rate of the polymers. Only the simplest mesoscopic polymer 
model, the Hookean dumbbell model plj . is discussed here. In the present study, the dumbbells do not react back on 
the shear flow. As it will turn out, turbulence is able to stop the stretching of the linear springs, but at a scale that 
is larger than the Kolmogorov scale of the flow.[27j Our approach intends to keep the polymer dynamics simple, but 
to take an advecting flow in its full turbulent complexity stemming from the Navier-Stokes equations. 

In chapter II we present the dynamical equations for the flow and the dumbbells and discuss characteristic param- 
eters. The dependence of local stretching rates on the shear rate is quantified in chapter III. The subsequent chapter 
summarizes our findings on the dumbbell dynamics where extension and angular statistics are discussed. Eventually, 
a summary and an outlook are given. 



II. MODEL AND EQUATIONS 
A. Advecting fluid 

The Navier-Stokes equations for a three-dimensional incompressible fluid are solved by a pseudo-spectral method 
using a second-order predictor-corrector scheme for advancement in time, pfij The equations of motion are 

-^ + (u-V)u = -Vp + ^V 2 u + f, (1) 

at 

V-u = (2) 

where u is the (total) velocity field, p the kinematic pressure field, and f the volume force density. The case without 
shear is solved in a three-dimensional box of side-length L = 2tt with periodic boundary conditions. The nearly 
homogeneous shear flow is modeled in a volume with free-slip boundary conditions in the shear direction y and 
periodic boundaries otherwise. Here the total velocity field follows by a Reynolds (de)composition as a linear mean 
part with the constant shear rate S and a turbulent fluctuating part 

u = (u) + v = Sye x + v . (3) 

The aspect ratio is L x : L y : L z = 2n : n : 2n. The applied volume forcing is a combination of an isotropic forcing 
that injects energy at a fixed rate of ej„ =0.1 and a shear forcing. Consequently, (e) = ej„ for the statistically 
stationary case at S = 0. More details on the volume forcing and the numerical scheme can be found in Ref. 29. In 
Tab. 1, we summarize some statistical properties of the flows that were studied. Since we want to resolve eventually 
a range of scales below 77, we are limited in the range of accessible Reynolds numbers although resolutions of up to 
512 x 257 x 512 grid points are used for the simulations. 



B. Hookean dumbbells 



The simplest description of polymer stretching can be accomplished by considering Hookean dumbbells. Their 
entropic elastic force is linearly dependent on the separation vector, R = x 2 — Xi, that is spanned between both beads 
at positions X2(£) and xi(i), respectively. When taking into account the elastic entropic force, hydrodynamic Stokes 
drag, and thermal noise the evolution equation for R reads 0, l30l | 



R = Au- 




where Au = u(x 2 , t) — u(x l5 i) is the relative fluid velocity at the bead centers. The last term is the thermal Gaussian 
noise with the following properties 



(&(*)> = 



(5) 
(6) 
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run 


l 


2 


3 


4 


5 


6 


S 





1/tt 


3/tt 


5/tt 


7/tt 


9/tt 


V 


1/30 


1/30 


1/30 


1/30 


1/30 


1/30 




0.100 


0.113 


0.285 


0.845 


2.268 


5.179 


R\ 


10.5 


9.1 


19.2 


29.8 


36.2 


43.7 


S* 





1.21 


4.85 


7.31 


7.58 


7.78 


St v 





0.17 


0.33 


0.32 


0.27 


0.23 


Tav/T 


20.5 


25.1 


18.9 


21.0 


112.8 


70.7 




8.39 


8.12 


6.44 


4.91 


3.84 


3.12 



TABLE I: Parameters of the direct numerical simulations. S is the constant mean shear rate, v the kinematic viscosity, (e) 
is the mean energy dissipation rate, R\ — Ak /5/(3(e}^)(i) 2 } is the Taylor microscale Reynolds number with the mean square 
of the turbulent velocity fluctuations (v 2 ), and S* = S{v 2 )/(e) is the dimensionless shear parameter. The averaging time T av 
is given in units of the large scale eddy turnover time T = (t> 2 )/(2(e)). The spectral resolution is indicated by k max r/ where 
k m ax = v2AT/3. All studies of the flow properties were done with N — 128. 



with i, j = x,y, z. It prevents the extension of a dumbbell to shrink below its equilibrium length 



with ks being the Boltzmann constant, T the temperature, and H the spring constant of the Hookean spring. Equation 
(|7|) follows from the equipartition theorem. Time constant r is the relaxation time of the polymers and is given for 
dumbbells as |2l| 

T = w> (8) 

where f = Qirpva is the Stokes drag coefficient with fluid mass density p and bead radius a. We sketch in Fig. 1 the 
coordinate system that is used. Since we restrict our study to Hookean dumbbells, the effects of finite extension are 
not imposed by the potential. Therefore the extended dumbbells experience both the smooth and rough scales of 
the advecting flow during their time evolution. Consequently, the velocity difference Au is not approximated by the 
linearization (R • V)u in Eq. Q as usually done. We solve the equations for the beads separately, but can reconstruct 
the full dynamics of R afterwards. An ensemble of up to 2 x 10 6 dumbbells, i.e., 4 x 10 6 beads, is advanced by a weak 
second-order predictor-corrector scheme simultaneously with the flow equations. 30J Initially, their center of mass is 
seeded randomly in space with a uniform distribution and an initial extension of Rq. All Lagrangian interpolations 
were done with a trilinear scheme since the velocity field is very well resolved as given by the spectral resolution factor 
kmaxV in Tab. 1. 



C. Stretching rates 

According to Eq. J2J the dynamics of dumbbells is subject to both, local stretching rate due to varying strain at 
small scales and restoring linear spring force. The separation vector between two fluid elements, |5r|, evolves as 

^Srj(t) = a jk (t) Sr k (t) for j,k = x,y,z, (9) 

where <Jjk(t) is the local stretching tensor along the Lagrangian trajectories. Consequently, 

5rj(t) = W jk (t,0)Sr k (0), (10) 

with the time-ordered exponential 

W jk (t, 0) = T+ exp Qf <T jk {r) dr^j , (11) 
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FIG. 1: Dumbbell coordinate system that is used throughout this text. It follows from the sketch that R x — RcostjjcosO, 
R y = i?sin (f/cos 8, and R z — sin 8 where R is the distance between both beads. The notation differs from conventional spherical 
coordinates, but has the advantage of giving alignment with the outer mean flow for <f> = 6 = 0. 4> is azimuthal angle and 8 the 
polar angle. 



where i5r(0) is the initial separation vector. The local stretching and contraction rates are quantified by the Lyapunov 
exponents 

1 \Sr^(t)\ 

A ' = A"«j <log( ]jT<4f>- (12) 

with i = 1,2,3 for three dimensions and Ai > A2 > A-j. |3lLl32T | The subscript L is an average over different realizations 
of the Lagrangian tracer tracks. The largest Lyapunov exponent, A1, quantifies the growth of the norm of the 
separation vector 5r. The exponents are long-time averages of Lagrangian simulations following the algorithm by 
Benettin et al. [3]]]. Figure 2 illustrates that the convergence of the Lyapunov spectrum for a larger shear rate 
requires a long time series. A sufficient convergence is reached there after about 70 large eddy turnover times the 
latter of which is defined as T — (v 2 )/(2(e)). Clearly this procedure has to be done for each S value anew. 

While the mean stretching time scale in a turbulent flow is given by the \± , the beads relax back to their coiled 
equilibrium distance within times of the order of r. The Weissenberg number which is defined as 

Wi = Ait , (13) 

is an appropriate measure for the competition of those mechanisms. For Wi < 1/2 the relaxation to equilibrium 
size is on average faster than stretching by the flow. Polymers are preferentially in a coiled state and their size 
distribution is stationary. In contrast, for Wi > 1/2 stretching proceeds in time until the finite extensibility limit 
is reached or turbulence stops the growth of dumbbells. This transition at Wi ~ 1/2 is known as the coil-stretch 
transition. l33| The pro perties of such a transition in random and short-in-time correlated flows have been studied 
analytically. |l7l ITsl Il9l |20| The probability density function (PDF) of end-to-end distance R = |R| exhibits then 
algebraic tails for Wi < 1/2, i.e. when R — + 00 one gets a scaling law P{R) ~ where q depends linearly on 

Wi -1 . In case of Hookean dumbbells, the PDF is not any longer normalizable at values Wi > 1/2. 

Although the analytical treatment is possible only for this limited class of random short-correlated flows, the coil- 
stretch transition appears in other turbulent flows. Experimental 0, Il5| and several numerical studies 0, Il3l Il9l 
l34l l35l l36l l3^ | all provide direct or indirect verifications of such a transition within different models for polymers 
or in distinct flow geometries. A sharp transition from a preferentially coiled to a preferentially stretched state at 
Wi = Wi c is however often not observable since different local flow topologies cause differently large mean extensions, 
as known from single molecule experiments. | l(l lll| 



III. STRETCHING PROPERTIES OF THE SHEAR FLOW 



Let us discuss now the stretching properties of the shear flow. The dependence of the Lyapunov spectrum on the 
shear rate S is shown in upper panel of Fig. 3. One can see that all three exponents grow monotonically in magnitude 
with increasing S. The lower panel of the same figure shows the largest Lyapunov exponent as a function of the Taylor 
microscale Reynolds number R\. Furthermore, the strength of the shear can be quantified by the dimensionless shear 
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FIG. 2: Time evolution of the three Lyapunov exponents Ai for the simulation with a shear rate S = 7/tt. The dashed line is 
the sum of the exponents which has to be zero due to incompressibility of the flow. Time is in units of the large scale eddy 
turnover time T = (w 2 )/(2(e)). 1000 Lagrangian tracer particles were seeded homogeneously at the beginning of the simulation. 
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FIG. 3: Upper panel: Lyapunov spectrum A* for i = 1, 2, 3 as a function of shear rate S. The inset replots Xi(S) in logarithmic 
axes and fits a power law Ai ~ S 3 ^ 2 to the larger shear rate values. Lower panel: Largest Lyapunov exponent Ai as a function 
of the Taylor microscale Reynolds number, R\ = ^5/(3(e)i/){v 2 ). 



parameter S* = S(v 2 )/(e). Table 1 indicates that this parameter saturates with growing shear rate to a value of 
about 8 while the Reynolds number is still increasing with growing shear rate. Such a saturation of S* is known from 
studies in homogenous shear flows in the stationary regime. |29L I38L l39| Our findings as given by Fig. 3 and Tab. 1 
underline the dependence of all turbulence parameters on the shear rate S. 

An analytic treatment has to rely on a simpler case. In Ref. 23, the fluctuating Kraichnan flow is added with a 
linear shear flow which gives a decomposition as in (J3J. The correlation matrix of the Gaussian fluctuating strain is 
defined as (ViVj(t)VkVi(t')) = DCijkiS(t — t'). It contains an amplitude parameter D that models the strength of the 
small-scale fluctuations and can be varied independently beside the shear rate S in the linear part. The amplitude D 
introduces a fluctuation time scale Tf — D^ 1 beside the shear time scale T s — S^ 1 . 
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In the large shear limit, which is defined as 

the molecules are most of the time aligned with the shear. The azimuthal angle 
(((>) ~ {D/S) 1 / 3, and satisfies the asymptotic power law distribution P{4>) ~ 
(</>} <C </> "C 1. The positive Lyapunov exponent follows to Ai ~ S(4>) and thus 

Ai ~ 5 2 / 3 . 

Relation 1|14|1 is known in Navier-Stokes turbulence as being close to the rapid distortion limit of a shear flow [icj . 
when the large shear rate determines the turbulent dynamics. 

There is a slight increase of Ai in comparison to the isotropic case for S < 1 as the inset of the Fig. 3 indicates. 
However, as the shear rate increases further we observe a crossover to a faster nonlinear increase with S. In fact 
Ai ~ y/ (e) jv and (e) ~ S 3 L 2 leads to the 3/2 scaling, 

c.3/2 



(14) 

t> (see Fig. 1) has a nonzero average 
l/4> 2 for an intermediate range of 

(15) 



Ai ~ S 6/z , (16) 

which differs from Q15[l . Although the range of our shear rates is limited, the observed scaling is consistent with this 
dimensional estimate. This argument and the saturation of the dimensionless shear parameter suggest that our shear 
flow operates in a different regime. Let us therefore compare shear and fluctuation time scale in the present case. 
The small-scale gradient fluctuations are captured by the Kolmogorov time t v = \J i//(e). T s has to be related to t v 
which is known as the Corrsin parameter. 41] Figure 4 shows that 



'/> r„ Jj^j or St v < 1 . 



(17) 



In fact l |lti|l dictates St v ~ S 1 " 1 / 2 . In order to support this behavior we add Tab. 2 where data from DNS at higher 
Reynolds number |29| and from turbulent boundary layer measurements [421 l43j are shown. The experimental data 
demonstrate that even close to the wall the Corrsin parameter does not exceed unity as can be seen for the data at 
y + = 34. We conclude that our studies of the polymer dynamics in the turbulent shear flow cannot be compared with 
the analytic results of Chertkov et al. |2^| 
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FIG. 4: Ratio of the fluctuation time scale to the shear time scale, St v . Present data are filled symbols and the open symbols 
are taken from Ref. 29 (see also Tab. 2). 



IV. POLYMER STRETCHING 



A. Statistics of dumbbell extension 



Let us turn now to the stretching of the dumbbells in the described turbulent shear flow. All parameters of the 
runs in Brownian dynamics are summarized in Tab. 3. The corresponding turbulence parameters are as in Tab. 1. 
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V 
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y + 


R\ 
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St v 


DNS1 


1/300 


1/ir 


- 


55 


4.2 


0.14 


DNS2 


1/300 


3/tt 


- 


175 


7.7 


0.16 


DNS3 


1/300 


8/tt 


- 


309 


9.8 


0.09 


HFI1 


Uair 


1576 


34 


151 


37.8 


0.97 


TTT7TO 


Vair 


on 


i ryin 

io4y 


I 


2.8 


0.14 


DNWl 


&air 


1447 


709 


1156 


50.3 


0.17 


DNW2 


Vair 


151 


7826 


1739 


20.7 0.05 


DNW3 


Vair 


78 


22753 


1299 


14.3 


0.04 



TABLE II: Further data on the ratio Sr^ taken from Ref. 29 (see series II in Tab. 1 there). In addition, data for two boundary 
layer measurements are listed. Measurements were done by Knobloch and Fernholz at the Hermann-Fottinger Institute (HFI) 
and the German-Dutch Windtunnel (DNW) (for more details on the experiment, see Ref. 42). The current values listed here 
are calculated from a data analysis by Jachens.l43|| The free-stream velocity for HFI is (Too = 10 m/s and for DNW (Too = 80 
m/s. For the experiments, the shear rate is given in 1/s and the kinematic viscosity of air is v air = 1.5 x 10~ 5 m 2 /s. 



s 


Wi 


TV 


N d 


v/Ro 


Ro/A 


run 2p 1/7T 


0.26 


512 


2 x 10 6 


5.5 


2.0 


run 3p 3/-7T 


0.36 


512 


2 x 10 6 


4.4 


2.0 


run 4p 5/-7T 


0.59 


128 


5 x 10 5 


1.7 


1.0 


run 4pl 5/-7T 


0.26 


128 


5 x 10 5 


3.4 


0.5 


run 4p2 5/ir 


0.26 


128 


5 x 10 5 


1.7 


1.0 



TABLE III: Parameters for the shear flow simulations with the dumbbell ensemble. Other parameters of the simulations are 
as in Table 1. Nd is the number of dumbbells and TV the grid points in streamwise direction which corresponds to a resolution 
of N x x N y x TV 2 = TV x (TV/2 + 1) x TV. We also list the ratio of the viscous Kolmogorov scale to the equilibrium extension of 
the dumbbells and the ratio of the equilibrium extension to the grid spacing. The grid spacing is A = 2-k/N . 



According to Eq. IjRjfl . effects of mean shear will enter the Weissenberg number as well. A series of three simulations 
(2p, 3p, 4p) is conducted for which r is kept fixed and the Weissenberg number changes in correspondance with \i(S) 
as shown in the upper panel of Fig. 3. In addition, runs 4pl and 4p2 are at Wi — 0.26 equal to that of run 2p, but 
at a larger shear rate of S = This series will demonstate in brief another possibility of changing the physical 

parameters, namely keeping S fixed and varying Wi via the relaxation time r. Furthermore, we will study the Rq 
dependence of the extension statistics. Our focus anyhow will be on the first series. 

In Fig. 5, we plot the instantaneous snapshot of all the dumbbells that are stretched beyond lOry for the shear 
rate S — 5/ir (run 4p2). A first glance on the plot indicates no preferential orientation of the polymers and a rather 
complicated pattern in which they are arranged. This is exactly our motivation for a more detailed study of the 
angular statistics which will follow at the end of this section. 

First, we analyze the statistics of the end-to-end norm R. Similar to isotropic flows we detect the occurance 
of the coil-stretch transition. The numerical computation of the probability density function (PDF) p(R, t) for 
Wi < Wi c — 1/2 and Wi > Wi c is shown in Fig. 6 and Fig. 7 respectively. Note that the definition of the probability 
density function is chosen such that J °° p(R,t)dR — 1. Runs 2p, 3p and 4p2 correspond to the preferentially coiled 
state of the dumbbells, i.e. Wi < 1/2. The PDFs are stationary after a relaxation time that increases with Wi and 
display fatter tails with increasing shear rate S. The dependence on the equilibrium extension Rq is shown in the 
inset of Fig. 6. We find that the PDF remains unchanged, except the very far tails. For both cases, Ro was taken 
well inside the regular and smooth viscous range of scales of the flow which makes our finding plausible. 

Run 4p corresponds with the preferentially stretched case for Wi > Wi c . We plot a sequence of PDFs each taken at 
the times indicated in the legend of Fig. 7. Although the small scale stretching is on average dominating the Hookean 
spring force the PDF p(R) reaches a stationary state. One may loosely think of a mean field picture that the relative 
velocity between the beads scales as Au(i?) ~ XiR. Therefore under the condition that X±R > i?/(2r), i.e. when the 
Weissenberg number is larger 1/2, the average inter-bead distance should grow with time. 

As mentioned in the introduction, we determine the relative velocity between the beads by taking the full increment 
between the individual bead velocities. Consequently, the dumbbells get stretched ultimately to scales R > r) for 
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FIG. 5: (color online) Instantaneous snapshot of the dumbbell distribution in the turbulent shear flow at S = 5/n and Wi = 0.26 
(run 4p2). Only the dumbbells with R > 10r; are plotted. The scale IO77 is the sidelength of the square on the bottom of the 
frontside of the box. 

Wi > 1/2. At those scales the relative velocity between the beads scales as Au(i?) ~ R a with < a < 1. Hence 
the growth of separation stops as the stretching term becomes subdominant in comparison to linear relaxation by 
the spring force. We add two external scales of our system into the figure. The integral scale of turbulence L int (see 
Ref. 29 for definition) marks the average extension of the largest structures in the shear flow and therefore the end of 
the inertial range. The box size L x is also shown. It can be observed that the sudden drop in the slope of p(R) for 
later times takes place at R < L X) but at R > Lint- Note that the turbulent motion is completely decorrelated across 
scales R > Li nt and that large extensions (which go even beyond R = L x ) become more and more unprobable. 

Such a scenario is not new, it is a main ingredient of the cascade model of drag reduction in bulk turbulence 
suggested by Tabor and De Gennes.|27| For a well established cascade with Au(R) ~ ((e)R) 1 ^ 3 a saturation scale of 
i?„ ~ Wi 3 ^ 2 ^ would follow. The mean extension (R) for the later plots of Fig. 7 is found to be at about 3rj, i.e. the 
saturation takes place on average beyond the viscous scale of turbulence. In our case, the saturation is a result of the 
nonlinear advection term in Eq. {3J in combination with the Hookean spring force. 

B. Anisotropy of stretching 

It is important to observe that the effect of shear cannot be merely translated to an effective scaling of the Weis- 
senberg number. The shear flow is not isotropic and hence the stretching shall not be so. To track the impact of shear 
we have plotted in Fig. 8 the PDFs of the individual components Ri for two cases below Wi — 1/2. As expected, the 
stretching in the streamwise direction becomes dominant when the Weissenberg number increases. The extension in 
y direction is always bounded to L y due to the free slip boundary conditions. This effect is also existing for very low 
shear rates. 

Figure 9 shows the second order moments (RiRj) for run 3p at Wi = 0.36. Effects of anisotropy are well detectable 
in the figure. The moment (R x R y ) has a finite positive value while the other two mixed moments were found to vary 
around zero. Moreover, the moments (R 2 ,) , (R y ) and (R 2 Z ) are no longer collapsing as for a lower shear rate. The 
stretching in the streamwise direction gives not only the largest extensions, but also the largest fluctuations. When 
comparing it with a graph of the turbulent kinetic energy of the corresponding run 3p (see inset of Fig. 9), we recognize 
similar patterns of time variation. They indicate that the large scale variations of the turbulence fluctuations have an 
effect on the stretching history, most dominantly in streamwise x-direction. For that reason we have to take statistical 
averages over longer time intervals with increasing S. 

Such fluctuations are known for homogeneous and nearly homogeneous shear flows. We recall that the scale of 
variation of the mean velocity gradient is outside the considered simulation domain. Simulations which start with 
isotropic initial conditions will show a self-similar growth of integral scale and turbulent kinetic energy until length 
scales of the finite simulation domain are reached. A statistically stationary regime is then accompanied by fluctuations 
of the turbulent kinetic energy. [2^, [3^, |3t| Such velocity (and kinetic energy) fluctuations show up as streamwise 
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FIG. 6: Probability density function (PDF) of the dumbbell extension R as a function of the shear rate S. The length R = |R| 
is rescaled by the corresponding equilibrium extension Ro of the advecting shear flow. Symbols are indicated in the legend. 
The slopes fitted to the tails are -3.0 for run 2p at Wi = 0.26 (solid line) and -1.5 for run 3p at Wi = 0.36 (dashed line). Run 
4p2 at 5* = 5/tt and Wi = 0.26 is also shown in the main figure (see Tab. 3). Inset: Ro dependence of the PDF of dumbbell 
extension for two different equilibrium length values as indicated in the legend (see also Tab. 3). 




FIG. 7: (color online) PDF of the dumbbell extension for the stretching at S — 5/tt and Wi = 0.59 (run 4p). The extension R 
is rescaled again by the corresponding equilibrium extension Ro. Curves at several times are plotted as indicated in the legend. 
The inclined arrow to the left indicates growth in time. T = (u 2 )/(2(e)) is the large scale eddy turnover time. The vertical 
dotted lines are for the integral scale of turbulence Lint and the box size L x — 2ir. 



streaks in the shear flow. Our studies are also in agreement with recent works by Stone and Graham 13] and Terrapon 
et aL0] who found preferential stretching near streamwise streaks. Both cases considered a so-called minimal flow 
unit that captures important features of a turbulent buffer layer. In the latter work such preferential stretching study 
was detected by an eigenvalue analysis of the local velocity gradient along the dumbbell tracks. The streaks decay 
due to linear instabilities and become less coherent when the Reynolds number increases. 0] We expect that such 
pronounced variations in the streamwise extensions will decrease for larger Reynolds numbers. 
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FIG. 8: Anisotropy of the stretching as quantified by the PDFs p(Ri/Ro) for the streamwise component i = x (upper panel), 
the shear component i — y (mid panel), and the spanwise component i = z (lower panel). 



C. Angular statistics 



The results of the last chapter suggest an investigation of the statistics of the azimuthal angle <j) which measures 
the orientation of the dumbbells with respect to the mean flow component. The azimuthal angle is calculated as (see 
also Fig. 1) 

4> = arctan (j^J , (18) 

with values between — n/2 and ir/2. Figures 10 and 11 show our analysis for the runs 2p, 3p and 4p (see also Tab. 3). 
To resolve the anisotropy of stretching, the contour plot of the joint PDFs of extension and azimuthal angle, p(R, (f>), 

is numerically computed and plotted in Fig. 10. They are normalized such that J °° dR f^j^ d<fi pC^>0) = 1- The 
lowest shear rate at S = 1/tt shows practically no variation with respect to 4> (see upper panel of Fig. 10). For highly 
stretched dumbbells the (f> probability distribution is essentially flat for almost all <p values exept a narrow region 
around (f> ~ 0. However this tendency of alignment with the mean flow is very weak and only visible for the most 
extended polymers. The dumbbell dynamics is very similar to that of an isotropic flow. 

The picture changes in case of higher shear rates as can be seen in the lower panel of Fig. 10. We observe a larger 
number of dumbbells aligned around the mean flow axis at <p — even when they are not very extended. Furthermore, 
in comparison with lower shear, it is less probable to observe dumbbells perpendicular to x-axis, i.e. at (f> = ±n/2. 
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FIG. 9: Time evolution of the polymer conformation tensor components as quantified by {RiRk} / Ro- Data are for S = 3/ir and 
Wi = 0.36. The symbols of the four components are given in the legend. Times are rescaled by the large scale eddy turnover 
time T = (v 2 )/(2(e)). Inset: Turbulent kinetic energy as a function of dimensionless time t/T. The long time analysis was 
conducted at N = 128. 



A closer inspection of both plots unravels an asymmetry in <p probability distribution. Dumbbells with <p £ [ — tt/2, 0] 
are unstable with respect to the mean flow in contrast to those with <p £ [0, 7r/2]. 
The PDF of the azimuthal angle is calculated by integration of the joint PDF 

P {4>) = / dRp(R,<f>). (19) 
Jo 

The angular distribution is shown in the upper panel of Fig. 11. The asymmetry between both quadrants is quantified 
by the following measure 

A{<f>) =p{4>) -p(—</>) , (20) 

with (f> G [0, 7r/2]. The measure A{<j>) is plotted in the lower panel of Fig. 11. Note that the integral of A{<j>) determines 
the total asymmetry of the PDFs. We find that with increasing shear rate the asymmetry of the angular distribution 
grows. The plot o{p(<f>) shows a growing peak around <fr = with increasing Wi. Fluctuations in the vicinity of <\> — 
are enhanced while the tails for very large 4> are depleted. 

The second angular degree of freedom is the polar angle defined as, 

9 = arctan f Rz \ , (21) 

where 9 G [— 7r/2,7r/2]. The statistics of 9 is therefore reconstructed by having the full information of the three 
components R^. The simultaneous statistical information of end-to-end norm R and the polar angle 9 is expressed by 
the joint probability distribution p(R, 9) which is normalized as J °° dR f^j^ ^ 0) = 1. 

Figure 12 shows the contour plots of the joint distribution p(R,6) for runs 2p and 3p. Now, 9 G [0, 7r/2] because 
the dynamics is invariant under the transformation 9 — > — 9 (which we verified in our simulations). The probability 
distribution of R at a fixed polar angle is given by vertical cuts in the two-parameter plane. For S = 1/n (upper 
panel of Fig. 12), the contour levels in the vertical direction are very close and insensitive to the polar angle. A low 
Weissenberg number at this shear rate is consistent with close spacing of the contours in R direction and indicates a 
rapid decrease of the extension probability. Moreover, the decay of p{R 1 9) with R at fixed 9 is more or less similar 
for the majority of polar directions. Recall that a larger polar angle 9 means that dumbbells are less aligned with 
the x-y shear plane. Consequently, orientations aligned with the shear plane and off that plane are equally probable. 
For S — 3/ 7T (lower panel of Fig. 12), the contour level spacing in the R direction is larger, obviously a result of the 
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FIG. 10: Joint probability density function p(R/Ro, <f>) for two different simulations. The contour levels of p(R/Ro, 4>) are given 
in units of the decadic logarithm and decrease in steps of 0.5. Upper panel: data are for run 2p with a shear rate S = 1/n and 
Wi = 0.26. Lower panel: data are for run 3p with a shear rate S = and Wi = 0.36. The contour lines for 10~ 3 and 10~ 2 
are indicated by exponents -3 and -2, respectively. 



increase of Weissenberg number. The other feature is that the projected distributions along the R direction decrease 
more rapidly with increasing the polar angle 9. This manifests the increase of the anisotopy in 9 direction. Therefore 
the dumbbells aligned with x-y shear plane are highly stretched only. 

To demonstrate this effect, an instantaneous snapshot of the most extended dumbbells for run 4p2 is shown in 
Fig. 13. Again, dumbbells with R > 10rj are plotted only, but now in two different sectors of the polar angle. A 
preferential alignment of the dumbbells with smaller polar angles - sector < 9 < 37r/20 is compared with sector 
77r/20 < 9 < 7r/2 - can be clearly detected from both panels (see also Fig. 5). 

Figure 14 shows the probability distribution of polar angle, 

/>oo 

p(9) = / dRp(R,6). (22) 
Jo 

The comparison of run 4p with runs 2p and 3p shows that the number of dumbbells aligned with the x-y plane 
increases with shear rate. This is reflected in the increasing maximum of the probability distribution around 9 w 0. 



V. SUMMARY AND DISCUSSION 



We present numerical studies of the stretching of Hookean dumbbells in turbulent flows with different mean shear 
rates which are described by the Navier-Stokes dynamics. The flow at hand is the simplest statistically stationary 
turbulent shear flow, a nearly homogeneous shear flow with a linear mean profile. It can be thought as the flow in the 
bulk of a turbulent shear flow with no-slip boundaries, e.g. in a turbulent plane Couette flow. Beside the turbulence 
structures above the Kolmogorov scale rj, the dumbbell dynamics around their equilibrium distance Rq < rj is well 
resolved. The latter is important for the lower Weissenberg numbers where the majority of dumbbells remain below 
the viscous scale. This constraint clearly limits the accessible range of Reynolds numbers. 

In the first part, we studied the stretching properties of the flow that are quantified by the largest finite-time 
Lyapunov exponent. For larger values of S our results infer that Ai ~ S* 3 / 2 which is consistent with a simple 
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FIG. 11: Upper panel: probability density function p(<j>) for the three runs 2p, 3p, 4p as given in Tab. 3. Lower panel: 
Corresponding asymmetry of p(<f>) as given by relation 1201 . 



dimensional analysis. Our studies and further supplementary data from other shear flows show that the Corrsin 
parameter does not get much larger than unity (in fact our data are all below 1). This point is important since it 
determines which flow time scale has to be compared with the relaxation time of the polymer chains. In our study 
the relevant time scale is ~ t, ( , in contrast to S 1-1 as in a recent analytical model, 

In the second part of the work, we studied the shear rate dependence of the stretching properties of Hookean 
dumbbells. The growth of anisotropy with shear rate is quantified by measuring the statistics of R4, the azimuthal 
angle tp, and the polar angle 9. When S grows a preferential orientation at azimuthal angles that are slightly larger 
than <f> = is observed. Moreover the probability distribution of the polar angle shows a preference for smaller values 
of 9 with increasing the shear rate. Therefore more polymers align with the shear plane as the shear rate gets higher. 
When keeping the shear rate fixed, the alignment of the dumbbells with the shear direction depends on the magnitude 
of the relaxation time r. For higher Weissenberg numbers and therefore larger relaxation times a more pronounced 
anisotropic stretching is then observed. Furthermore, the p(R) is found to be insensitive with respect to a variation 
of the equilibrium length Rq for most extension scales. 

Streamwise streaks which are a characteristic feature of shear flows are found to cause variations in the polymer 
extension fluctuations. In agreement with previous studies for the case of a minimal flow unit of a turbulent boundary 
laver [9l . these structures are found to enhance the stretching. For Wi > 1/2 the Hookean potential the inter-bead 
separations extend to the spatial scales where they experience a rough relative velocity difference. Therefore we detect 
a saturation of the stretching, even for Weissenberg numbers beyond 1/2 and in agreement with cascade models for 
drag reduction in the bulk of turbulence. [27j 

It will be interesting to investigate the behavior of finitely extensible nonlinearly elastic (FENE) dumbbells in our 
flow configuration. A further parameter, the maximum length Lq of the dumbbells, will have then an impact on the 
results. In that case well resolved velocity data below the Kolmogorov scale become very important. This is due to 
the fact that even for the largest polymer chains in turbulent flows Lq < rj. These studies are part of our future work. 



14 



50 



40 




20 40 60 80 




9 in deg 

FIG. 12: Joint probability density function p(R/Ro, 9) for two different simulations. Data and line spacing are the same as for 
Fig. 10. Exponents for levels 10 -3 and 10 -2 are indicated again. 
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